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We model a Superconducting Single-Electron Transistor operating by repulsive interactions. The 
device consists of a ring of Hubbard clusters, placed between electrodes and capacitively coupled 
to a gate potential. In each cluster, a pair of electrons at appropriate filling feels a weak efi'ective 
interaction which leads to pairing in part of the parameter space. Thus, the system can host many 
bound pairs, with correlation induced binding. When the charging energy exceeds the pairing energy, 
single-electron tunneling prevails; in the opposite regime, we predict the Coulomb blockade pattern 
of two-electron tunneling. This suggests that in tunneling experiments repulsion-induced pairs may 
behave in a similar way as phonon-induced ones. 

PACS numbers: 73.23.-b, 73.40.Gk, 81.07.Nb 



I. INTRODUCTION 

In recent years, a variety of transport experiments have 
been reported in molecular size systems, such as quan- 
tum dots and nanotubes, as a contribution to the current 
boost towards of the progress in nanoscale technology. 
From the theoretical side, circuits of several kinds have 
been modeled jj, and in applied electronics the Single- 
Electron- Transistors are among the most important 
devices. These are realized by connecting a nanoscopic 
conducting island to metallic leads and to a gate voltage. 
The energy gaps existing between states with different 
number of particles allow to fix the number of electrons in 
the island very sharply; as a consequence single electrons 
can tunnel to or from the conductor. Even more appeal- 
ing situations arise when the above scenario is compli- 
cated by electron-electron interactions, as in the case of 
a Superconducting-Single Electron Transistor (S-SET). 

A S-SET is a mesoscopic device obtained by linking ca- 
pacitively a superconducting grain to two normal leads 
and to a gate electrode as The latter allows one 

to control the number N of electrons on the grain by 
tuning the gate voltage Vg. Such a system has been 
studied both experimentally0| and theoreticallyQ|^0 
in great detail during the past years. In a normal island 
the parity of N oscillates between even and odd values, 
by varying Vg] conversely in a superconducting island N 
is always even because of the paired nature of the ground 
state. Therefore the S-SET transport properties in the 
linear regime are governed by Andreev reflection under 
the critical temperature Tc of the central island, while 
above Tc single electron tunneling prevails. This leads 
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to well pronounced Coulomb blockade peaks of the con- 
ductance G = dI/dV\v=o as a function of the gate volt- 
age. In particular the parity-controlled tunneling pro- 
duces 2e/Cg periodic peaks in the pair-tunneling regime, 
in contrast with the e/Cg periodicity of the normal sys- 
tem (here Cg is the capacity of the gate electrode). This 
behavior is well reproduced by models 's'i'g'i'?! using a 
gate controlled BCS Hamiltonian Hbcs', the connection 
to free electron leads employs a tunneling Hamiltonian, 
usually treated by second-order perturbation theory. 

In the present article we propose a model for a S-SET 
with a strongly correlated, repulsive Hubbard-like model 
instead of Hbcs as the "superconducting" grain Hamil- 
tonian. That is, we look for a superconducting response 
entirely driven by the electronic correlations rather than 
by the phonon-mediated effective attraction. The occur- 
rence of two-electron tunneling in non BCS systems was 
observed by Ashoori et a/.Q in a l/xm GaAs tunnel ca- 
pacitor. Purely electronic mechanisms were proposed to 
explain this behavior and the GaAs quantum dot models 
ranged from a semiclassical description^ to a Hubbard 
model framework (l^. Unlike the systems considered by 
Refs.0 |lOl |. in our gedankenexperiment , like in a S-SET, 
the tunneling current is due to many bound pairs hosted 
by the device in a wide range of gate potentials. 

The plane of the paper is the following. In the next 
Section we introduce the microscopic model that we are 
going to study. Section Hill is devoted to determine some 
important properties of the strongly correlated central is- 
land. We show that the electronic correlations provide a 
non-trivial characteristic energy which can be compared 
with the electrostatic charging energy in order to dis- 
tinguish between a normal regime and a superconduct- 
ing one. In particular in these two regimes the parity 
of the number of particles in the ground state oscillates 
exactly like in a S-SET. In Section Hvl we explicitly cal- 
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FIG. 1: (a): Scheme of the strongly correlated S-SET. The 
device consists of A Hubbard clusters arranged in a ring and 
linked symmetrically to one another, (fe): Pictorial represen- 
tation of Ha and Ht . 



culate the conductance as a function of the gate voltage 
by using a master equation approach 'll]. It is found 
that the linear response of our strongly correlated device 
shows Coulomb blockade pattern. A normal behavior is 
observed in the non-correlated and in the very strongly 
correlated regimes; while in the intermediate case, the 
spacing between the conductance peaks doubles. Finally 
the conclusions are drawn in Section FVl 



II. THE MODEL 



Let us consider the grand-canonical Hamiltonian 



Ti — fiNtot — -ffdcvicc + H] 



leads 



(1) 



Here i?dcvicc is an extended Hubbard model of the central 
island coupled capacitively to a gate voltage Vg] i?ioads 
describes the left and right reservoirs, supposed to be 
identical free electron gases for simplicity; Ht is the tun- 
neling Hamiltonian, that connects the central device to 
the leads; /i is the chemical potential and Ntot is the total 
number of particles operator. The carriers are electrons, 
of charge — e, e > 0. Let us examine these three terms in 
detail. As in previous workfl^. the central island consists 
in a ring of A identical 5-site centered-square Hubbard 
clusters, denoted by the index a, see Figs^, b. Each 
cluster is described by the Hamiltonian 



4 



with the creation operators on the a-th cluster p^^ for 

the central site, and p]j j^, i = 1, .., 4 for the remaining 4 

sites; fia^ia — pjj i^Pa.ia-, and cr is a spin index. 

In the device, each cluster a is linked to the two near- 
est neighbors ones (denoted by a -I- 1 and a — 1) by the 
hopping Hamiltonian Hr (see Fig^) whereby a particle 
in the i-th site of the a-th cluster can hop towards the 



i-th site of the /? = a ± 1-th clusters: 

A 4 

Hr = rJ2 E E^^'Up/^.^ + I^-^-)- (3) 

a—l (3—a±l a,i—l 

^^dovicc also contains an electrostatic charging energy 
term due to an effective capacitance C of the central is- 
land. Finally the island is connected capacitively to the 
gate which is at a potential Vg (see Fig^). Therefore we 
have 

Q — 1 

where N is total number of particles operator in the cen- 
tral device. We remark that the capacitive term is essen- 
tially long-ranged and accounts for the monopole contri- 
bution to the charging energy, while the U terms depend 
on the way the charges are distributed in the island. In 
all electrostatic terms N should be referred to an average 
population corresponding to a neutral situation; but, ac- 
tually, any shift N N — (N) would produce a constant 
and a linear term in N that just modifies fi. 

Both leads are free electron gases with chemical poten- 
tials /i^, 7 = l,r; hence 



with fi — fii — fir ~ ey|l6| where V is the bias. 
Finally the tunneling Hamiltonian is taken to be 



(5) 



r^^)(ct,,/[A/2],,.+h.C.) 



(6) 



where the are eigen-operators of the noninter- 

acting term of Ha- Y^ta^iP^afiaPa^Ky + pi.iaP"''^'^'l = 

Y.n.a^v fi,-na fa,^^ Wc obscrvc that the tunnel junc- 
tions connect two opposite clusters to the leads; namely 
the a = 1 cluster is linked to the left electrode and the 
a = [A/2] cluster to the right lead (here [x] means the 
integer part of x), see Fig^. Note also that T'-'''' is inde- 
pendent of ry, in other terms we are using "white" wires, 
that is, leads that do not filter electrons according to 
the square symmetry of electronic states in each cluster. 
This is a simple way to ensure the 3D nature of the leads, 
which is essential to allow Andreev refiection. 

In the next Section we draw some relevant properties 
of i/dovicc which mimic the behavior of Hbcs despite the 
presence of strong electronic correlations. 



III. PROPERTIES OF THE CENTRAL ISLAND 

In order to understand the physics of the device that 
we propose, it is useful to foucus first on properties of 
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FIG. 2: A (in eV) as a function of U/t. The maximum 
binding occurs at (7 ~ 5t where A ^ -0.042 t. For U > S4t 
A becomes positive and pairing disappears. 



Ha, referred to a single 5-site cluster. The Hamiltonian 
Ha is a prototype example of electronic pairing from re- 
pulsion; this is signaled by the property A < 0, where 
A = £(4) + e(2) — 2e(3) and e{m) is the ground state 
energy with m electrons. There is pairing aX m — A for 
U /t < 34, the minimum value is A = — 50meV at t = leV 
and U ~ 5eV and the binding energy is |A| (see FigHJ. 
The mechanism has been investigated elsewhere|13i| ,14] . 
and need not concern us here; we just say that broadly 
speaking it is a lattice counterpart of the Kohn-Luttinger 
mechanismjlSl] . 

2 

When a negative A occurs, its competition with ^ 
determines the parity of the number Ngs of electrons in 
the ground state of -ffdcvicc for small values of the inter- 
cluster hopping T. 

Here we are interested in the behavior of the device 
at low temperatures kT <^ |A|. Up to o(t) and o(t^) 
corrections, the ground state energy E{N) of the central 
device with fixed even or odd number of particles TV is 



Ae(2) -I- (iV - 2A)/p 



•■ 2 



2C 



A)A- 



for even TV 



En = <! 



Ae(2) + (iV - 2A)/p + (f - A - i)A 



^Nef 
2C 



e{Vg - ^i)N 



for odd N 



(7) 

where Ip = e(3) — £(2). Since one bound pair exists at 
m = 4 electrons in the 5-site cluster, the first bound pair 
in the A-cluster system appears at = 2A-I-2 electrons. 
From Eq.(|7I), it follows that in the range 2A < A^ < 4A, 
Ngs is always even if the pair binding energy overcomes 
the charging energy. We call the situation when |A| > 

2 

|t the superconducting regime. Otherwise the system is 
normal and any A'^ is lowest in a range of V^. In FigOwe 
plot the ground state energy of i?device as a function of 
Vg in both regimes, for r = 0. It also relevant to focus on 
the critical values of Vg where ground states of different 
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FIG. 3: iJjv — Es, where is the ground state energy of 
i^dcvicc, versus Vg, for several numbers N of particles, (a) in 
the superconducting regime; (b) in the in the normal regime. 
Solid lines are used for even N and dotted lines for odd N. 
In the insets we plot the grand-canonical averages of N in 
Hdovicc We used A = i, t = leV, r = 0, kT = O.OOleV, 
C = 50e/V; with this choice e^/C = 0.02eV. In (a) U = 5eV 
(A = -0.043eV ) and in (&) U = 0.2eV (A = -O.OOOSeV). Vg 
is in V, En is in eV. 



Ngs cross. We define (AVg)„ the spacing between the 



critical values in the normal regime and {AVg)s 
spacing in the superconducting regime; it holds 

(A^,)„ = # + f 

{AVg)sC - t . 



such a 



(8) 



The charge fluctuations in the superconducting regime as 
a function of the gate voltage are about double spaced 
with respect to the normal regime, that is the typical 
condition experimentally realized in a S-SEtQ. 

We can also visualize the previous results by plotting 
the gran-canonical average of the number of particles in 
the isolated central device as function of Vn . We use the 



standard definition 



z 



(9) 



where Z = Tr[e-(^<i— -^*)/'=^]; Tr is dominated by the 
low energy states of TJdevicc with r = 0. In Fig^ one 
can observe the so called Cooper staircase, characteristic 
of the superconducting regime. 

Below, for computational convenience, we assume r ^ 
|A| and deal with Ht perturbativelvfl^ . So, the critical 
values of Vg where level crossing occur, are spread into 
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FIG. 4: (A'^) versus Vg.((i) in the superconducting regime; (b) 
in the in the normal regime. The parameters are the same as 
in FigOl Vg is in V. 



intervals of width o(t) in the normal regime and o(r^) 
in the superconducting regime. Anyway the qualitative 
behavior is still close to Figs l3l4l 



IV. CALCULATION OF THE LINEAR 
CONDUCTANCE 



Next, we consider the effects of small bias voltage 
V applied between leads, i.e. the linear conductance 
G = dl/dV for F — 0, versus the gate voltage Vg. In 
the present article, we follow the approach proposed by 
BeenakkerfTlj and get the formula for the conductance 
from a master equation. We take T^'^T^'"' << kT « 
A I in order to provide that (i) the parity of the ground 
state is stable with respect to thermal effects, (m) the 
elementary tunnel processes between the leads and the 
central devices involve few particles at a time and the 
broadening of the levels of iJdevice due to the presence 
of the leads is smaller than the thermal one. As dis- 
cussed by Beenakkerjllj. these limitations characterize 
the Coulomb blockade regime. 

In the normal regime single-electron tunneling domi- 
nates. The theory works very much like in Ref.[ll| and 
we calculate the first-order transmission rates 



p(7) ^ ^7^(7)^ 



E 

?7,(j 



{(lN-l\foL.,ri(7\^N) 



(10) 



from the iAr_i-th state of the central device with N par- 
ticles (denoted |«Ar)) to the qjv_i-th state of the central 
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FIG. 5: Linear conductance G'^^' as a function of Vg in the 
pair-tunnehng regime. We used A = 4, t = leV, U = 5eV, 
T = 0.0005, kT = O.OOleV, C = 50e/V, ^ =0. Vg is in V. 



device with A'^ — 1 particles ( \qN-i)) via tunneling to 
the left (1) or rigth (r) lead. Remember that a = 1 for 
7 = I and a = [A/2] for j = r. In the actual calculations 
below, we obtain [zat) in second-order perturbation the- 
ory in Ht- ; we mix the degenerate ground-state multiplets 
of energy En of iJdcvice, which determine the low-energy 
properties of the system. By first-order perturbation the- 
ory in Ht, one gets the familiar formula: 



4A 

pe_ 

KT ^ ^ 

7V=2A+1 j2v,9iv-i 

xp(o)(zw)[i -/(£;,„ 



(11) 



Here, P^^^{iiy) is the Boltzmann equilibrium probabil- 
ity for occupying the eigenstate |iiv) with energy 
/ is the Fermi distribution function, and p is the den- 
sity of states at the Fermi level in both leads. Each 
term in Ea. lfTT)l depends on Vg through the statistical fac- 
tor p(°)(zAr)[l — f{Ei^ - Eqj^_-^)] ar id p roduces the well 
known Coulomb blockade behavior [llj . The linear 
conductance is highly suppressed unless the gate voltage 
is fine tuned at Ei^ ~ ^'gw-i ~ Mi where sharp peaks of 
occur. The second-order contribution in Ht depends 
on o([r''']^), rates which are negligible with respect to 
the odT^]"^) r coefficients; therefore we can safely avoid 
working out the second-order current in this regime. 

When |A| > |t, only even N have an important weight 
in the appropriate range of Vg (see Fig|31a); therefore the 
resonance condition Ei^^ = Ej^ -^ — p, never holds and the 
first-order conductance G'-^-' is highly suppressed for any 
value of V^. In this pair tunneling regime, accordingly, 
we must go on calculating the conductance up to second- 
order in iyT0@[3- Three-body, four-body transitions 
and so on can be disregarded, however, as T^'^ and T'-'"' 
are both small compared to the charging energy. Since 
electrons can get paired in the device but not in the leads, 
we may think of the second-order processes in terms of 
Andreev reflections. First one of the two electrons tun- 
nels from one lead to the device (which is in the \mN-2) 
state) and forms a virtual excited state. Then the second 



5 



12 


A=0 U=0 




1 n 






^ 6 




/y 11 
c 1 1 




■ (a) 


Ac 


4 




2 



^^^N=13 J\^N=14 J 


\n=15J \ 



0.27 0.28 0.29 0.3 0.31 0.32 0.33 0.34 



17.5 
15 
12.5 
10 
7.5 
5 
2.5 




A=0 U/t=34 

r 

c 

(b) " 




0.4 0.425 0.45 0.475 0.5 0.525 0.55 



FIG. 6: Linear conductance G*-^' as a function of Vg in the 
single-electron-tunneling regime: (a) U = 0, {b) U = 34eV. 
The other parameters are the same as in Fig|3 Vg is in V. 



peaks for every Vg such that Ei^ = Era^_2 +2/x, while the 
conductance is strongly suppressed elsewhere. Each peak 
has a correlation weight due to the coefficient Si^^mj^ ^j 
containing all the microscopic information on the corre- 
lated ground states of the central device. 

For illustration, we numerically computed the conduc- 
tance for a central device with A = 4. In the supercon- 
ducting regime, as a function of the gate voltage 
is shown in FigElfor t = leV, U = 5eV, r = O.OOOSeV, 
kT = O.OOleV, C = 50 e/. Note that r < | A| = 0.043eV, 
and |A| > e^/C = 0.02eV. G^"' shows neat peaks, with 
spacing (AVg)sc = 2e/C = 0.04 V. 

This superconductor-like behavior depends on the exis- 
tence of pairing. As a countercheck, we calculate the lin- 
ear conductance G^-* in the normal regime, when A = 0. 
We can obtain this condition in two ways, namely in the 
non-interacting case when [/ ~ and in the very strongly 
correlated regime, when U ~ 34eV, the other parameters 
remaining the same as in FigEl G(^) is plotted in FigEl 
since we are mainly interested in the period of the reso- 
nances, we use constant F's and plot the results in arbi- 
trary units. Indeed, for A = the period of the resonant 
peaks is e/G ~ 0.02V, i.e. a half of the period in the 
superconducting case. 



one tunnels into the device and form a bound pair (|iAr) 
state). 

In principle in second-order, one should also take into 
account the cotunneling processes |l8j. which leave the 
population of the central island unchanged. Such pro- 
cesses provide a current away from the resonances. Any- 
way, as long as our device is in small bias and low tem- 
perature regime, the cotunneling current is found to be 
negligible, as in the case of Ref..6j. Therefore the se- 
quential tunneling — > Jat is the major transport 
mechanism, and it is possible only at the two-electron 
degeneracy points. The full derivation of the solution of 
the detailed balance equations will be presented else- 
where. Setting 



(II) _ 1287reV' t('^'t(-)' 



r'y'-'-i — 



one gets for the conductance 
G(" 



(12) 



4A,cvcn 



(11) 



^ ^ ^ ^ ^iN ,mN -2 

N=2A+2 iN,mN-2 



g 2kT 



sinh[(S,„ - E„,^_, - 2n)/2kT] 



(13) 



The amplitude 



X! X! (^N\fi^,^\lN-l){lN-l\fl^^mN-2) 
l„-l a,p.,ri.v 



takes into account the second-order process governing the 
Andreev reflection. Ea. (|13|) predicts Coulomb blockade 



V. SUMMARY AND CONCLUSIONS 

We have shown that the Hamiltonian Ti, models a S- 
SET in the linear regime. We pointed out that the 
repulsion-induced pairing occurring in ifdovice fixes a 
characteristic energy |A| which competes with the elec- 
trostatic charging energy e^/G. As in any S-SET, there is 
a normal regime, where |A| < e^/G and a superconduct- 
ing regime, where |A| > e^/G. In the flrst case the parity 
of the electron number in the ground state oscillates be- 
tween even and odd values and the transport properties 
are governed by single-electron tunneling. Conversely in 
the superconducting regime the parity is always even and 
the major transport mechanism is sequential tunneling of 
pairs. The explicit calculations have been performed for 
a ring of four 5-site clusters, but a general expression for 
the linear conductance is also derived. 

Our results suggest a systematic way to produce a well 
controlled periodic two-electron pattern, even without 
any conventional superconductivity' an array of quan- 
tum dots similar to the one in Ref. |31 could be designed 
to this purpose. 

Finally we underline that the model we propose is very 
flexible with respect [i) to the size and the shape of the 
Hubbard clusters, (m) to the topology of the cluster ar- 
ray forming the central device. Indeed a wide variety 
of Hubbard clusters show the A < property at proper 
fillings, which is actually the key feature at the basis of 
our device; we could construct many alternative devices, 
based on graphs with different topologies, also in view of 
possible single-electronics applications to more complex 
circuits than a transistor. 
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